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ABSTRACT 

We study the outer density profiles of dark matter haloes predicted by a generalized 
secondary infall model and observed in a dissipationless cosmological simulation of 
a low-density flat cold dark matter model with the cosmological constant. We find 
substantial systematic variations in shapes and concentrations of the halo profiles 
as well as a strong correlation of the profiles with the environment in which the 
haloes are embedded. In the iV-body simulation, the average outer slope of the density 
profiles, (3 {p (X r~^), of isolated haloes is /3 « 2.9, and 68% of these haloes have 
values of (3 between 2.5 and 3.8. Haloes in dense environments of clusters are more 
concentrated and exhibit a broad distribution of f3 with an average value higher than 
the average (3 for isolated haloes. For haloes located within half the virial radius of 
the cluster from the center values /3 « 4 are very common. Contrary to what one may 
expect, the haloes contained within groups and galaxy systems are less concentrated 
and have flatter outer density profiles than the isolated haloes: the distribution of /3 
peaks at « 2.3 — 2.7. The slope f3 weakly anticorrelates with the halo mass Mh- The 
concentration decreases with Mh, but its scatter is roughly equal to the whole variation 
of this parameter in the galaxy halo mass range. The mass and circular velocity of 
the haloes are strongly correlated Mh oc V^ with a « 3.3 and « 3.5 for the isolated 
haloes and haloes in clusters, respectively. For Mh « 10^^/i~^Mq the rms deviations 
from these relations are AlogM^ — 0.12 and 0.18, respectively. Approximately 30% of 
the haloes are contained within larger haloes or have massive companions within three 
virial radii. The companions are allowed to have masses larger than ~ 0.3 the mass of 
the current halo. The remaining 70% of the haloes are isolated objects. We find that 
the distribution of f3 as well as the concentration-mass and Mh — Vm relations for the 
isolated haloes agree very well with the predictions of our seminumerical approach 
which is based on a generalization of the secondary infall model and on the extended 
Press-Schechter formalism. 
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1 INTRODUCTION 

The dark matter (DM) haloes are thought to be objects 
within which luminous galaxies form and evolve. Thus, the 
properties and evolutionary features of the observed galaxies 
should be related to their haloes. According to the hierar- 
chical scenario, the DM haloes form via collapse of primor- 
dial density fluctuations in the expanding Universe. Cosmo- 
logical A'^-body simulations provide a direct way to study 
this process. Nevertheless, only recently the simulations be- 
came accurate enough to resolve the internal structure of the 
galaxy-size haloes produced in these simulations. Navarro, 



Frenk, & White (1996, 1997; hereafter NFW) found that 
for a large range of masses and background cosmologies the 
density profiles of equilibrium haloes are well approximated 
by the universal profile: 



p{r) 



r/rs{l -\-r/raY 



(1) 



where r^ is a characteristic radius, and ps is a character- 
istic density. NFW found that these two parameters are 
connected, and that the remaining free parameter depends 
only on mass. This free parameter is the concentration 
cnfw = rv/ra, where r„ is the halo virial radius. This con- 
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elusion was later extensively tested and confirmed with other 
numerical simulations (e.g., Cole & Lacey 1996; Kravtsov, 
Klypin & Khokhlov 1997; Moore et al. 1998; Bullock et al. 
1999). Studies based on analytical and seminumerical meth- 
ods also have shown that eq. (1) is a good approximation for 
the density profiles of typical equilibrium cold dark matter 
haloes (e.g., Syer & White 1998; Avila-Reese, Firmani, & 
Hernandez 1998; Salvador-Sole, Manrique & Solanes 1998; 
Raig, Gonzalez-Casado & Salvador-Sole 1998; Henriksen & 
Widrow 1999; Nusser & Sheth 1998; Kull 1999; Lokas 1999). 
Nevertheless, the applicability of the NFW profile has some 
limits. Avila-Reese et al. (1998) (hereafter AFH98) have 
shown that the NFW profile describes well only the structure 
of typical haloes formed from the most probable hierarchi- 
cal mass aggregations histories (MAHs). Depending on the 
MAH, diverse density profiles are possible. There is another 
limit on the applicability of the NFW profile: it is formally 
valid only for isolated haloes, which are in equilibrium and 
which do not have large companions. These limits on the 
applicability of the NFW profile point out to the necessity 
to find in the numerical simulations what kind of profiles 
have the haloes that deviate from the NFW case and how 
the different environments in which the haloes are embedded 
influence on their density profiles. 

Recently, Jing (1998,1999) has studied the density pro- 
files of hundreds of DM haloes in high-resolution A^-body 
simulations. He concluded that although the NFW profile 
describes the structure of a large fraction of haloes in equi- 
librium, there are other haloes whose density profiles deviate 
from the NFW shape. In particular, the haloes with signif- 
icant internal substructure show large deviations from the 
NFW profile. Bullock et al. (1999; see also Jing 1998, 1999) 
have found a substantial scatter in the parameter cnfw that 
may be due to the fact that the NFW profile does not de- 
scribe well some of the profiles. Bullock et al. also find that 
Cnfw depends on the environment. 

To determine the shape of the most inner parts of the 
halo density profiles the simulation should have a very high 
resolution. Whether this shape is r~^ as for the NFW profile 
or not, it is still matter of debate (Kravtsov et al. 1998; 
Moore et al. 1998,1999; Jing 1999). The resolution of the 
numerical simulations discussed herein is not sufficient to 
address the issue of the central slope. 

In this paper we study intermediate and outer regions 
{r^Vs) of the halo density profiles of a large sample of DM 
haloes obtained in a cosmological simulation of a flat low- 
density cold dark matter model with cosmological constant 
(ACDM). Details of the simulation are described in Colin et 
al. (1999) and Kravtsov & Klypin (1999). To identify haloes, 
we use the Bound Density Maxima (BDM) halo flnding al- 
gorithm (Klypin et al. 1999). The main difference with the 
results presented by Jing (1998, 1999) is that we study both 
the isolated haloes and the satellite galaxy-size haloes re- 
siding inside groups and clusters. As the result, we are able 
to examine the environmental distribution of the galaxy-size 
DM haloes and the influence of the environment on DM halo 
structural properties. 

Analytical and seminumerical methods, alternative and 
complementary to the expensive cosmological N-ho&y sim- 
ulations, have proved to be useful. These approaches allow 
one to follow some particularly chosen features and phenom- 
ena of DM halo formation. Here, the results obtained from 



the A'^-body simulations will be compared and, when possi- 
ble, will be interpreted in the light of the predictions of the 
seminumerical method presented in AFH98. This method is 
based on a generalization of the secondary infall model (see, 
e.g., Zaroubi & Hoffman 1993) and uses the extended Press- 
Schechter formalism in order to calculate the hierarchical 
MAHs. 

In §2.1 we briefly describe the numerical simulation and 
the algorithm used to identify DM haloes and to obtain their 
density proflles. An outline of the seminumerical approach is 
given in §2.2. The spatial distribution of the identifled haloes 
according to the environment is presented in §3. In Section 4 
we examine the outer density profiles of the haloes and the 
way in which these profiles change with the environment. 
The structural correlations of the haloes as a function of the 
environment are analyzed in section 5. The concentration- 
mass and mass- velocity relations are presented in §5.1 and 
§5.2, respectively. The estimated dispersions for these rela- 
tions and the correlation among them are also presented. 
In § 6 we discuss some of the results. The summary and 
conclusions of the paper are given in §7. 



2 NUMERICAL AND SEMINUMERICAL 
SIMULATIONS 

2.1 Af-body simulations and halo identification 
algorithm 

Structure of DM haloes only slightly depends on cosmol- 
ogy (e.g., NFW; Cole & Lacey 1996; Kravtsov et al. 1997; 
Avila-Reese 1998; Firmani & Avila-Reese 1999a). Therefore, 
the results for a representative cosmological model should 
be sufficient to outline the general behavior and trends of 
the structural properties of the DM haloes. In this study, 
we use a flat cold dark matter model with the cosmologi- 
cal constant (ACDM). The model has the following param- 
eters: the density of matter is fio = 0.3, the density due 
to the vacuum energy is Oa = 0.7, the Hubble constant 
is Ho = lOO/i km s^^Mpc"^ with h = 0.7, the amplitude 
of perturbations on Sfe^^Mpc scale is as = 1. The numer- 
ical simulation was done using the ART code described in 
Kravtsov et al. (1997). Details of the simulation are pre- 
sented in Colin et al. (1999). The simulation followed evo- 
lution of 256"^ particles in a 60/i^^Mpc box. The mass of a 
DM particle is 1.1 x 10^/i^^Mq. The peak force resolution 
is lAh~^ kpc. High mass and force resolution are very im- 
portant for survival of DM haloes in dense environments of 
groups and clusters of galaxies. 

The Bound Density Maxima (BDM) halo identifica- 
tion algorithm (Klypin et al. 1999) was applied to find the 
DM haloes. The algorithm locates maxima of density within 
spheres of radius 10/i~^ kpc and then removes unbound par- 
ticles. The algorithm produces a catalog of DM haloes con- 
taining coordinates, velocities, and density profile of bound 
particles for each halo. The density profile is used to find 
the maximum circular velocity Vm = {GM/rY^"^ , radius and 
mass of the halo. We define the halo radius r^ as the min- 
imum of the virial radius r„ and the truncation radius rt- 
The former is defined as the radius at which the average 
density of the system is Ac{z) times the average density of 
the universe at redshift z, where /S.c(z) is determined from 
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the spherical collapse model. For the model we use here 
Ac{z = 0) = 334 (e.g., Bryan & Norman 1998; but see a 
recent paper by Shapiro, Iliev, & Raga 1999 where a more 
proper and self-consistent treatment of the spherical collapse 
was carried out; for an Einstein-de Sitter universe they ob- 
tained that Ac{z = 0) is « 11% smaller than the standard 
value) . The truncation radius is the radius where the spheri- 
cally averaged outer density profile flattens or even increases. 
This radius marks the transition from the halo to the sur- 
rounding environment. Only a small fraction of the haloes 
identified in our simulation (--^ 6%) have rt < r„. The frac- 
tion of truncated haloes is larger for the non-isolated haloes 
than for the isolated ones. The mass of the DM haloes Mh 
is defined as the mass enclosed within rn- 

The BDM algorithm is capable of finding haloes with 
20-25 bound particles. In the simulation there were 9073 
identified haloes with this lower limit on number of parti- 
cles. Analysis of these haloes indicates that haloes with max- 
imum circular velocity Kn larger than 90 km s~^ (6819) are 
not affected by numerical effects and/or details of the halo 
identification (Gottlober, Klypin, & Kravtsov 1999). Nev- 
ertheless, because we need to find the shape of the density 
profile, we restrict ourselves to more massive haloes with 
Vm > 130 km s^'^. The number of haloes in our final cata- 
log is 3498, which still is large enough for our purposes. All 
haloes in the catalog have more than 200 particles. 

2.2 The seminumerical method 

AFH98 and Avila-Reese (1998) presented an approach to 
study the gravitational collapse and virialization of DM 
haloes formed from the Gaussian density fluctuations. The 
first step of the method is to generate hierarchical MAHs 
of DM haloes. We use the extended Press-Schechter approx- 
imation based on the conditional probabilities for a Gaus- 
sian random field (Bower 1991; Bond et al. 1991; Lacey & 
Cole 1993). For a given present-day mass, we generate a 
set of MAHs using Monte Carlo realizations. We follow the 
aggregation history of the main progenitor by identifying 
the most massive subunit of the distribution at each time. 
Then, the gravitational collapse and virialization of the DM 
haloes formed with the MAHs is calculated assuming spher- 
ical symmetry and adiabatic invariance during the collapse 
with an iterative seminumerical method. 

This method is based on the secondary infall model 
(e.g., Zaroubi, & Hoffman 1993). This model is modified to 
allowing non-radial motions and arbitrary initial conditions 
(MAHs in our case). The only free parameter in this ap- 
proach is the ellipticity of the orbits eo — rperi/rapo, where 
rpori and rapo are the pericentric and apocentric radii of an 
orbit, respectively. The parameter eo mainly influences the 
central structure of a halo: the more circular are the orbits 
(larger eo), the shallower is the inner profile. A''-body sim- 
ulations indicate that eo is typically 0.1-0.4 in cluster-size 
haloes (Ghigna et al. 1998). Here, we set eo equal to 0.15. 
This is the value for which the density profile of a halo of 
10^^ Mq produced with an average MAH has the same pro- 
file of an isolated well resolved halo of the same mass found 
in our A'^-body simulation. 

In order to start the Monte Carlo realizations, we fix 
the present-day halo mass, referred here as the nominal 
mass Mnom. At any time, the outer shells that encompass 



this mass are still in the process of virialization. The mass 
shells that are already virialized roughly correspond to those 
within the virial radius r„ at which the mean overdensity 
drops below the critical value Ac{z) given by the spherical 
collapse model. Analyses of haloes identified in numerical 
simulations show that at radii smaller than r„ the mater 
is indeed close to a virial equilibrium (e.g., Cole & Lacey 
1996; Eke, Navarro, & Frenk 1998). At radu between r„ and 
2r„ the matter is still falling onto the halo, while at larger 
radii, the matter is expanding with the universe. The mass 
contained within r„ is the virial mass Mv which, depending 
upon the MAH, is equal to 0.7-0.9 times Mnom (see also KuU 
1999). Because in the numerical simulations the mass of the 
haloes is defined by M„ (only in a few cases is defined by 
the mass at the truncation radius), for the seminumerical 
simulation we also use M„. 



3 THE ENVIRONMENTS OF GALAXY DM 
HALOES 

We divide haloes into two broad categories. One category 
is constituted by haloes whose centers do not lie within 
the radius of any other halo of equal or larger maximum 
circular velocity. We shall call these haloes distinct. Note, 
however, that the distinct haloes may contain other smaller 
haloes within their radii. The other category of haloes is 
residing within radii of haloes of larger maximum circular 
velocities. The haloes of this category are further divided 
into three sub-categories according to the size of their par- 
ent halo. If the parent halo has maximum circular veloc- 
ity of Vm > 600 km s~\ 350 km s"^ < Kn < 600 km s"\ 
or Kn < 350 km s^^, we will refer to them as haloes in 
clusters, tn groups, and tn galaxies, respectively. The limits 
which define the circular velocities of the cluster, group, and 
galaxy haloes are arbitrary. Nevertheless, they refiect veloc- 
ity ranges of real clusters, groups, and galaxies. It should be 
taken into account that the max;imum circular velocity of 
galaxy-size systems typically increases by a factor of 1.2-1.4 
due to dissipation in the baryonic component (AFH98; Mo, 
Mao & White 1998). 

Distinct haloes may or may not have massive neighbors. 
We shall call isolated those haloes that do not have a large 
companion with Kn°™^ > fvVm within 2>rh, where Vm and 
rh are the maximum circular velocity and radius of the cur- 
rent halo, and we have fixed factor fv ~ 0.7. In §5.2 we find 
that halo mass is related to the circular velocity approxi- 
mately as Mh oc Vm^'^- The constraint on circular velocity 
of the companion corresponds roughly to M™™'' > 0.3Mh, 
where Mh is the total mass of the current halo. Thus, an iso- 
lated halo is an object not contained within other halo and 
without massive companions up to a relatively large dis- 
tance. If a halo is not contained inside other halo (distinct) 
but has at least one massive companion (K™°™'' > 0.7Vni) 
within 3rh, we consider it as belonging to a multiple system. 
In fact, most of the multiple systems (~ 80%) are just pairs. 
That is why we shall refer to this class as the haloes in pairs. 

Table 1 gives the numbers and percentages of galaxy 
satellite haloes, haloes in clusters, and haloes in groups. 
Only 12.5% of all haloes belong to the category of haloes 
contained inside larger haloes. This fraction remains almost 
constant if we include in our catalog smaller haloes with 
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Kn < 130 km s^^ and with less than 200 particles. With the 
aim to find the fraction of haloes not contained inside larger 
haloes (distinct) but with massive companions (multiple or 
pair systems) , we analyze the surroundings of each of the dis- 
tinct haloes in search for companions. We may ask ourselves 
what is the distance dj^^p to the nearest companion with 
y^omp y o.7Kn, where Kn is the maximum circular velocity 
of the current halo. In Figure 1 we present differential and 
cumulative distributions of dSmp normalized to the radius 
rh of the halo. Although we usually consider only haloes 
with Vni > 130 km s~^, companions were allowed to have 
smaller circular velocities (> 90 km s^^). This was done to 
allow even a small halo (Vm ~ 130 km s~^) to have a chance 
to have companions as small as 0.7 of their own circular ve- 
locity. We find that only a small fraction (~ 2%) of haloes 
of the category in study contains a halo of mass larger than 
--0.3 of their mass (K™™^ > 0.7Ki) within their total ra- 
dius Th- Most of the haloes of this category have companions 
with y™™^ > 0.7Kn as far as 2-4 times their radius. The 
isolated haloes, as defined above, constitute the 80% of the 
distinct haloes and the 70% of all the haloes (see Table 1). 

The large fraction of isolated haloes found in the nu- 
merical simulation actually strongly depends upon the pa- 
rameter fv- In Figure 2 we plot the fraction of isolated 
haloes with respect to all the haloes as a function of fv- 
Because we prefer to limit our catalog only to haloes with 
Kn> 90 km s~^ (this is the minimum velocity allowed for 
the companion haloes) , the limit on Vm of the isolated haloes 
has to be increased when fv decreases in order for the sam- 
ple to remain complete. That is why as fv decreases we 
should use catalogs with larger limits on Vm- The number 
of isolated haloes significantly decreases when the minimum 
mass of the companions decreases. Our results agree with 
the halo-halo correlation function for isolated haloes. 

The fractions of objects in different systems found at 
2: = in our numerical simulation roughly agree with what 
is observed in the Universe: 60% — 70% of galaxies are in 
the field (most of them are disc galaxies), 30% — 40% are 
in groups (e.g., Ramella, Geller, Huchra 1989; Nolthenius, 
Klypin, & Primack 1994), and 5% — 10% are in clusters 
(Bahcall 1988). It should be noted that some of the pair 
haloes (~ 17% of all haloes in the sample) might be classified 
as small groups composed of two relative large galaxies and 
a few small satellites. Moreover, as it was mentioned above, 
~ 20% of these haloes actually have more than one massive 
companion, i.e. they form multiple systems. 



4 DENSITY PROFILES 

4.1 TV-body simulations 

The mass resolution in our simulation (m = 1.1 x 
lO^/i~^M0) is not sufficient to resolve central parts of most 
of our haloes. As the result, we focus on the structure of the 
outer profile. Our first question is whether the halo density 
profiles have the shape of the NFW profile. AFH98 (see also 
Jing 1998,1999) find that the DM haloes actually have a 
range of density profiles for a given mass where the average 
profile may be described by the NFW profile. Using the large 
sample of haloes identified in our numerical simulation, we 
fit the spherically averaged density profile of the haloes by 
the following function: 



p(r) = 



-(1 + — )'3- 



(2) 



This is a generalized NFW profile where the slope /3 of the 
outer part of the profile (p(r) oc r~^ for r » rs) may be 
different from the slope /3 = 3 of the NFW profile. 

The density in the inner regions with less than 50 par- 
ticles has a rather large shot noise. Therefore, we use only 
those bins which have more than 50 particles inside them. 
For most of the DM haloes, the radius from which the halo 
has more than 50 particles is ~ 0.3 — 0.8 of the radius where 
the maximum circular velocity r^ is reached (for the NFW 
profile the radius rm is about 2.2ra). Since our interest is 
in galaxy-size haloes, we additionally restrict the sample to 
haloes with Vm < 350 km s~^. This reduces the number of 
haloes to 3347 (out of 3498). 

The frequency distribution of the parameter /3 obtained 
for the galaxy-size haloes is plotted in Figure 3. We have 
found that (3 does not depend on mass; there is actually an 
indication for a weak anticorrelation. The arrows in the hor- 
izontal axis, from left to right, indicate the 16%, 50% and 
84% of the cumulative distribution, respectively. In other 
words, roughly 68% of the DM haloes have values of /3 
between 2.50 and 3.88, where the median corresponds to 
/3 « 2.94. This result does not depend strongly on the qual- 
ity of the fit. In Figure 3 the frequency distribution of (3 for 
those profiles that were fitted with an accuracy better than 
(x^/^bins)^''^ < 7%0 is plotted with the thin solid curve. 
Approximately 33% of all the haloes used in this analysis 
have a fitting with this accuracy. The number of bins A'^bins 
varies from halo to halo depending on its size. This is why 
we divide x^ by Abins in order to have an estimator of the 
goodness of the fit. The distribution becomes only a little 
narrower and slightly shifts to smaller values of /3 than in 
the case when all the profiles are considered. 

An error analysis of the slope /3 is important in our 
case because the number of particles in an average halo is 
not very large. We roughly estimate limits of the error in the 
following way. For a given halo, we generate a set of density 
profiles drawn from an ensemble of profiles with the mean of 
the original halo profile and with deviations defined by the 
Poisson noise due to finite number of particles. The set is 
used to estimate the errors in l3 produced by finite number 
of particles in the halo. Specifically, for each radial bin of a 
halo profile we find the number of particles in the original 
halo and then perturb it assuming the Poissonian distribu- 
tion of the particles in the bin. Repeating this procedure for 
every bin and several times for each halo, we get a set of 
density profiles for a given halo. Applying the fitting pro- 
cedure to each of these density profiles, we obtain a set of 
values for the fitting parameters (/3, Vs, and ps) for which we 
can estimate the standard deviations. This method provides 
a way to estimate the uncertainty on /3. 

Having in mind that /3 does not depend on the mass, 
we have applied the experiment to three groups of haloes 
with the slope /3 around 2.50, 2.94 and 3.88. These values 
of the slope correspond to the 16%, 50% and 84% of the 



* In our case, the quantity x^ used for the minimization of the 
fitting is relative (dimensionless) because we fit the logarithm of 
the density ( x^ = J2i O'^sifti/Pan))'^, where pa„ and pi arc the 
analytical and measured values of the density, respectively 
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cumulative distribution of this parameter. For each group 
we selected dozens of haloes, and for each halo we applied 
the Monte Carlo experiment 30 times. The average standard 
deviations of the parameter /3 for each group are shown in 
Figure 3. Note that the dispersion of l3 increases with /3. This 
is expected because for haloes of a fixed mass, the number of 
particles in external bins is smaller when /3 is larger. There- 
fore, the Poisson noise of external bins is larger for larger /3. 
The dispersion of the other two parameters, particularly po, 
have the opposite trend. Thus, a significant contribution to 
the uncertainty in l3, particularly when l3 is large, is proba- 
bly due to the relatively small number of halo particles (see 
also Figure 4a) . A minor contribution to the uncertainty in /3 
may be due to the fitting technique. The Monte Carlo exper- 
iment that we have applied to the haloes, also can be viewed 
as a procedure to produce small deformations in the density 
profiles. In a few cases the fitting technique can give com- 
pletely different values of /3 for a set of these profiles. These 
cases typically happens when the density profile abruptly 
changes from a very shallow slope to a very steep slope. In 
these cases the scale radius rs is typically fixed at very large 
and unphysical values. 

4.2 Dependence of the outer halo density profiles 
on environment 

The outer part of the density profile eq. (2) is described 
by the parameter /3. We find that this parameter depends 
on the halo environment. In Figure 4 the distribution of 
P is shown for haloes in different environments. The fre- 
quency of haloes is defined with respect to the number of 
objects in the given category (environment). Because most 
of the haloes are isolated, the distribution for all the haloes 
presented in Figure 3 remains almost the same for the iso- 
lated haloes. The distribution corresponding to haloes with 
Vm > 350 km s^^ plotted in Figure 4a is slightly narrower 
and it has lower amplitude at high values of /3. But the dif- 
ference is small. Thus, it appears that large and galaxy-size 
haloes have similar distribution of the slope /3. The differ- 
ences are probably due to the fact that more massive haloes 
have more particles, and, thus the scatter on the outer den- 
sity profiles for them is smaller, particularly when /3 is large. 

The external slope for pair and galaxy satellite haloes 
in most cases is /3 ~ 2.2 — 2.6 (Figure 4b), which is shallower 
than the slope in the NFW profile. For the haloes in groups 
the distribution of /3 is wider and shifted to larger values of 
/3 with a maximum frequency around /3 = 2.6 — 2.8. In the 
case of galaxy haloes in clusters, the distribution of /3 is even 
wider than in the other environments, with a maximum in 
/3 = 3.1 — 3.4. If we select only the galaxy haloes contained 
within half the total radius of the clusters, then we find 
that values of /3 ~ 4.0 — 4.4 are more frequent (dashed line 
in Figure 4c). It should be considered, however, that the 
uncertainty in the determination of /3 is large when l3 is large 
(see the error bar that accounts for the average standard 
deviation estimated for the cluster haloes with values of l3 
near 3.9). In any case, the trend of the parameter /3 with 
the environment is clear. 

With the aim to visually judge the quality of the fitting, 
in Figure 5 we plot the spherically averaged density profiles 
and the corresponding fitting to eq. (2) for isolated haloes, 
galaxy satellite haloes, group haloes, and cluster haloes. In 



this Figure we have plotted for each category three randomly 
chosen haloes, each one from a given range of masses, and 
with values of /3 around the corresponding maximum of its 
distribution. Except for a few cases, the fitted density profile 
describes very well the structure of the DM haloes. To show 
the quality of the fitting even in the cases where the uncer- 
tainties of the fit are high (when (3 is large) , we present the 
profiles for the haloes in clusters with /3 around 4.0 instead 
of around 3.3. In Figure 5, comparing the profiles of the most 
massive haloes (upper curves) with the less massive (lower 
curves), it can be appreciated how the number of particles 
influence on the quality of the result. The profiles of the less 
massive haloes (less particles) are noisier than those of the 
more massive haloes. Therefore, the fitting for the former is 
more uncertain than for the latter. 

We also use another way to fit the halo density profiles 
in order to check our results. Instead of leaving the param- 
eter l3 free (see eq. (2)), it was fixed to two different values 
/3 = 3 and /3 = 4, where the former corresponds to the NFW 
profile and the latter to the Hernquist profile (Hernquist 
1991). As already was mentioned, the quantity x'^Z-^bins may 
serve to some degree as a criterion of the goodness of the 
fit. In Figure 6 the values of (x'^Z-^bins)^ obtained for the 
NFW and Hernquist profiles are compared for two samples: 
the isolated haloes (panel a) and the cluster haloes (panel 
b). For the former, the NFW density profile in most of the 
cases is a better approximation than the Hernquinst profile. 
For the latter sample, the density profiles of a large fraction 
of haloes are better fitted by the Hernquist profile than by 
the NFW. This supports the result that haloes in clusters 
tend to have steeper outer density profiles than the NFW 
shape. 

4.3 Results from the seminumerical method 

Using the seminumerical approach we produce catalogs of 
halo profiles for each chosen mass Mnom. In order to esti- 
mate the slope l3 and the concentration we apply the same 
fitting procedure used for the results of the numerical sim- 
ulation. Because in the seminumerical approach we are not 
able to introduce dynamical effects related to the environ- 
ment like the tidal stripping, and due to the assumption of 
spherical symmetry, several effects related to non-sphericity, 
particularly the major mergers, are not considered. This is 
why the haloes produced in the seminumerical simulations 
correspond to the isolated haloes identified in the numerical 
simulations. In Figure 4(a), the dashed line represents the 
distribution of the parameter /3 obtained from the seminu- 
merical approach. As in the numerical simulations, in this 
case we also find only a very weak dependence of /3 with the 
mass. The Figure 4a shows that both the simulation and the 
seminumerical approach produce similar results. 

Strictly speaking, the distributions of /3 for the haloes 
produced in the numerical and seminumerical simulations 
presented in the upper panel of Figure 4, do not correspond 
to the same estimate. In the case of the numerical simu- 
lations, as it was discussed in §4.1, the uncertainty due to 
the relatively small number of particles introduces an ex- 
tra scatter on the distribution of /3. The procedure we have 
used to calculate this scatter likely overestimates the errors 
because it does not preserve mass of the halo. Moreover, in 
few cases some contribution to the calculated scatter can be 
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also due to the ambiguity of the fitting technique (see §4.f ). 
This also applies for the haloes in the seminumerical simu- 
lation. The intrinsic distribution of /3 for the isolated haloes 
obtained in the numerical simulation will be narrower than 
that presented in the upper panel of Figure 4. Nevertheless, 
we expect that it does not differ much from the distribution 
obtained in the seminumerical simulations. 

According to the seminumerical approach the differ- 
ences in the structure of the haloes are mainly related to 
the dispersion of the MAHs. The origin of this dispersion is 
due to the statistical nature of the primordial density field. 
Haloes that have larger rates of mass aggregation at earlier 
times (early collapse), have density profiles more concen- 
trated and typically have steeper outer slopes than haloes 
with larger mass aggregation rates at later epochs. We find 
that the outer slope /3 is particularly sensitive to the behav- 
ior of the MAH at late epochs (close to z = 0). For example, 
if the halo suffers a very late major merger, its outer profile 
slope P will be small (/3< 2.6). On the contrary, if the mass 
aggregation rate is very small for 2< 1, then (5 tend to be 
larger than 3. 



5 CONCENTRATION AND STRUCTURAL 
CORRELATIONS OF THE HALOES 

5.1 Mass vs. concentration 

For the CDM-like power spectra of fluctuations, the hierar- 
chical MAHs of the DM haloes are such that on average the 
less massive objects attain a given fraction of its present-day 
meiss slightly earlier than the more massive ones. Therefore, 
on average, less massive haloes are more concentrated. For 
the NFW profile, cnfw = r^/rs is a reasonable and physi- 
cally motivated parameter of concentration. However, in the 
case of the more general profile given by eq. (2), the scale 
radius Ts has different physical meanings for different values 
of /3. That is why it is desirable to define a concentration 
parameter independent from the fitting applied to their den- 
sity profiledTI. From the numerical and seminumerical simu- 
lations we find that the ratio between the halo radius and 
the radius containing 1/5 of the total mass. 



Cl/5 



'■1/5 



Th 



r(Mi,/5), 



(3) 



is a reasonable estimator of the halo concentration for most 
of the haloes. This is because, for the typical haloes (those 
with /3 ~ 2.7 — 3.0), ri/5 is near to the radius r^ where com- 
monly the mass profiles differ more from one to another. The 
parameter Ci/5 correlates with Cnfw for a given value of /3. 
In Figure 7 we plot these two parameters for halo profiles 
from the numerical simulations with /3 ~ 2.5 (stars), ~ 3.0 
(dots), and ~ 4.0 (crosses). This plot shows the limitation of 

T As a matter of fact, it is not possible to characterize the struc- 
ture of the DM haloes emerged from a stochastic density fluctu- 
ation field with only one (universal) parameter. The diversity of 
density or mass profiles associated with the diversity of MAHs 
(see e.g., AFH98; Firmani & Avila-Reese 1999a, b), certainly re- 
quires for their description more than one parameter. Neverthe- 
less, a good level of approximation may be attained with a min- 
imum number of parameters when the parameters are appropri- 
ately defined 



the parameter cnfw when the NFW shape is generalized to 
the profile given by eq. (2). Our results show that ci/5 and 
j3 are weakly correlated. For practical purposes Ci/5 and (3 
may be considered as two independent parameters; each one 
is associated with different characteristics of the mass distri- 
bution of the haloes. While (3 describes the density profiles 
at the outer regions, Ci/5 deals with the overall mass distri- 
bution down to intermediate (« 0.5 — l.Or^) radii. 

In Figure 8 we plot the parameter Ci/5 as a function of 
virial mass for the isolated haloes (a), the haloes in groups 
and galaxies (b), the haloes in cluster (c), and the haloes 
obtained in the seminumerical simulations (d) (in this case 
Th = r„ always). The average values of Ci/5 at each mass bin 
were used in these plots. We find that the dispersions of C1/5 
have an approximate normal distribution. The standard de- 
viations corresponding to each bin are presented in Figure 
8 with the dashed lines. The dispersion in the concentration 
is related to the dispersion in Vm- According to the mass- 
velocity relation (see §5.2), a halo with ~ 200 particles (our 
lower limit) have Vm ~ 125 — 130 km s~^. Nevertheless, due 
to the dispersion, there are haloes of this mass with smaller 
or larger circular velocities. Thus, in order to avoid statisti- 
cal incompleteness in the estimate of the dispersion of the 
concentration at masses closer to the lower limit, we have 
fixed the lower limit on velocity to 100 km s^^ instead of 
130 km s~^ . The thin solid lines in each panel are the linear 
regressions to all the haloes of the corresponding sample. We 
find that the low mass haloes tend to be more concentrated 
than the high mass haloes. The haloes in clusters, although 
with a significant dispersion, also tend to be more concen- 
trated than isolated haloes. This is in qualitative agreement 
with Bullock et al. (1999). Since Ci/5 is a parameter com- 
pletely independent of the fitting, the fact that haloes in 
clusters have larger values of Ci/5 than the isolated haloes of 
the same mass, suggests that the former have steeper density 
profiles than the latter. 

Figure 8 shows that the standard deviations of the con- 
centration Ci/5 are of the order of the whole variation of 
this parameter with mass in the galaxy-mass range. This 
variation should be taken into account by analytical and 
semianalytical works on galaxy formation and evolution. 

The Ci/5 — Mh relation predicted by the seminumeri- 
cal approach is in excellent agreement with the results for 
the isolated haloes. Nevertheless, we should note that this 
agreement is expected. As was pointed out in §2.2, in the 
seminumerical approach we have to fix a parameter, eo, re- 
lated to the ellipticity of the particle orbit. Here, we have 
fixed eo to the value for which the density profile of a model 
of Mh ~ W^^ Mq (produced with the average MAH) agrees 
with the profile of a typical isolated halo of the same mass. 
Note, however, that the trend of Ci/5 with the mass and its 
scatter predicted with the seminumerical approach, are inde- 
pendent from the normalization. If we fix eo using the profile 
calculated with the procedure outlined in NFW (1997) then 
the concentrations Ci/5 are smaller than those obtained with 
the profiles of our numerical simulation by a factor ~ 1.3. 
Indeed when we fit the density profiles from the numerical 
simulation to the NFW profile (eq. (1)) and define the virial 
radius in the same way NFW did it (the radius where the 
average density of the halo is 200 times the critical den- 
sity), we obtain that the average values of the parameter 
Cnfw are ~ 1.2 — 1.4 larger than the values calculated with 
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the NFW 1997 procedure. Recently, using high-resolution 
Ai'-body simulations, Moore et al. (1999) have also reported 
values of cnfw 50% higher than the values given in NFW97. 



5.2 The mass-velocity relation and its dispersion 

The average Vm corresponding to several mass bins for iso- 
lated haloes (a), cluster haloes (b), and the haloes obtained 
in the seminumerical simulations (c), are plotted vs. the 
mass in Figure 9. The dashed lines represent the respective 
standard deviations, and the thin solid lines are the linear 
regressions for all the haloes of the corresponding sample. 
As previous works have shown (e.g., NFW97; AFH98; Bul- 
lock et al. 1999), a strong correlation of the form Mh oc V^ 
at galaxy scales is found. The average slope a we find for 
isolated and cluster haloes is ~ 3.29 and ~ 3.50, respec- 
tively. In the case of the seminumerical simulations (to be 
compared with the isolated haloes), a ~ 3.22. Again, the nu- 
merical and seminumerical approaches give a similar result. 
The Mh — Vm relation exhibits a dispersion which is due to 
the statistical nature of the primordial density fluctuation 
field (AFH98; Avila- Reese 1998). 

In Figure 10 we present the fractional rms deviations 
of the velocity ay/ < Vm > as a function of mass for the 
isolated haloes and the haloes in clusters as well as for the 
haloes in the seminumerical simulations. The haloes in clus- 
ters have larger deviations than the isolated haloes. Note, 
however, that due to the small number of haloes in clus- 
ters, the noise in the determination of their deviations is 
high. In the clusters the haloes are subject to tidal stripping 
which is able to change the original structural properties of 
the haloes (Klypin et al. 1999). These changes introduce an 
extra scatter in the mass-velocity relation. The deviations 
obtained with the seminumerical approach are very simi- 
lar to those of the isolated haloes. The results for a ACDM 
model of Eisenstein & Loeb (1996), who used a very sim- 
plified method to estimate the deviations of the M^ — Vm 
relation, are also very similar to those obtained here from 
the numerical simulations. The fractional rms scatter in ve- 
locity may be translated into the logarithmic rms deviation 
of the mass: AlogM^ = alog(l -I- ay/ < Vm >), where a is 
the slope of the Mh - Vm relation. For Mh ^ lO^^/i'^M© we 
obtain AlogM^ =0.12, 0.18, and 0.11 for the isolated and 
cluster haloes, and for the haloes from the seminumerical 
simulations, respectively. 

In Figure 11 we present the correlation among the resid- 
uals of the Mh — Vm and C1/5 — Mh relations for the isolated 
and cluster haloes. For a given mass, the more concentrated 
are the haloes the larger are their Vm- The scatter in this 
correlation is larger for the haloes in clusters than for the 
isolated haloes. Some haloes deviate from the correlation 
among the residuals; they apparently have too large Vm 
for their concentrations. In fact most of these haloes are 
those which were truncated (§2.1); their masses and radii 
are smaller than the virial mass and radius, while their ve- 
locities and radii where is contained 1/5 of the mass remain 
nearly the same. As one may see from Figure 11 the loga- 
rithmic deviations from the Mh — Vm relation are roughly a 
factor 2 smaller than the corresponding deviations from the 
Ci/5 — Mh relation. 



6 DISCUSSION 

6.1 The density profiles of haloes in clusters and 
groups 

In §§ 4.2 and 5.1 it was shown that the outer shape and the 
concentration of the halo density profiles appear to be infiu- 
enced by the environment. As in previous studies (Ghigna et 
al. 1998; Okamoto & Habe 1998; Klypin et al. 1999), we also 
found that the haloes in clusters typically have steeper outer 
slopes than the NFW profile. Naively, one could expect that 
haloes in groups have outer profile slopes fiatter than those 
ot the haloes in clusters, but still steeper than the slopes of 
the isolated haloes. Our analysis shows that this is not the 
case. The satellite haloes in groups and galaxy-size systems 
as well as the pair haloes, have typically fiatter outer profile 
slopes and are less concentrated than the isolated haloes. 
Therefore, the halo density profiles do not follow a continu- 
ous trend along the cluster-group-field sequence. This result 
suggests that the differences between clusters and groups 
can not be viewed as a simple sequence in density. 

Why does the outer density structure of the galaxy-size 
DM haloes depend on environment? Tidal stripping plays an 
important role for haloes inside clusters: haloes that have 
been subject to tidal stripping have steeper outer density 
profiles than the NFW profile, while the haloes recently ac- 
creted onto the cluster have profiles in agreement with the 
NFW profile (Ghigna et al. 1998; Okamoto & Habe 1998; 
Klypin et al. 1999). This might be the case for some of our 
haloes. The MAH also infiuences the structure of the halo. 
We find that many of the haloes in clusters could have outer 
density profiles steeper than /3 ~ 3 because they have more 
concentrated profiles than the isolated haloes, and this might 
be because they formed earlier than the latter. For example, 
intherangeof masses of 4x10" /i"^ M0-5xlO"/i"^ M©, 
the density in the central bins of haloes in clusters is typi- 
cally 1.5-2.0 times larger than in the case of isolated haloes. 
Since in both cases the mass is roughly the same, then the 
external profile slope should be steeper for the cluster haloes 
than for the isolated ones. It should be taken into account, 
however, that if the haloes in clusters were tidally stripped, 
then their original masses have been decreased by the strip- 
ping. This could also explain why the central densities of 
present-day haloes in clusters are larger than those of the 
isolated haloes of the same mass. It seems that haloes in 
clusters tend to have steep outer density profiles due to 
both effects: (i) because they formed earlier than the iso- 
lated haloes in such a way their density profiles result more 
concentrated than the profiles of the isolated haloes, and (ii) 
because their outer parts were affected by the tidal stripping. 
Note that in the latter case the original halo concentration 
Ci/5 should be smaller. This is because the total mass of the 
halo Mh decreases only as ~ Inr^ (roughly M{r) ex In r at 
the outer halo parts). Thus, the halo radius rh is truncated 
due to the tidal stripping while ri/5 remains approximately 
the same. 

In groups, which are smaller and less dense than clus- 
ters, the tidal stripping is not a significant process, and the 
typical epochs of formation of haloes in groups do not dif- 
fer much from those of the isolated haloes. Therefore, the 
structure of the haloes in groups is not substantially affected 
neither by tidal stripping nor by the epoch of formation. 
Probably, the effects of recent aggregation and interactions 
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between the group members are more important than the 
stripping for group haloes. For galaxy-size haloes in clus- 
ters the roles of the two processes are reversed (Okamoto & 
Habe 1998). The profiles of some haloes in groups could be 
shallower than the equilibrium NFW shape because the halo 
is caught just when it begins to share the particles with a 
nearby companion. It is also possible that, even after a long 
time of virialization, due to the merging of substructures, 
the particle orbits are more circular than in the case of "un- 
perturbed" haloes (isolated). Therefore, these particles do 
not penetrate to the central regions and the density profile 
is shallower than for the unperturbed haloes. These situa- 
tions are even more probable for the galaxy satellite and pair 
haloes. The correct answer to the question of how and why 
do the structures of the DM haloes depend on the environ- 
ment have to come from a careful analysis of their evolution 
in different environments. This work is currently in progress. 



6.2 The origin of the mass-velocity relation 

The DM haloes exhibit a tight power-law relation M„ — Vm° 
between their masses and maximum circular velocities (§5.2) 
with the slope a « 3. It appears that the shape of the power 
spectrum of primordial perturbations is responsible for the 
slope. The power spectrum of fluctuations of the CDM mod- 
els is such that the concentration of the DM haloes only 
slightly depends on mass. Let us analyze the NFW profile 
(eq. 1), which describes well the density profiles of a large 
fraction of haloes in the numerical and seminumerical sim- 
ulations. The maximum circular velocity V-m of a halo is 
equal to: 



2 _ GM{< r„ 



(4) 



where 



2.16- 



cnfw 



(5) 



is the radius at the maximum circular velocity Vm ~ V{rm 
Integrating the NFW profile up to the radius r™ we find: 



M(< r„) ^ 0.467M„//(cNFw), 



(6) 



where f{x) = ln(H-a;) — a;/(l + a;), and the concentration pa- 
rameter cnfw ~ rv/rs is a weak function of the virial mass 
M„ (e.g., NFW). According to the definitions introduced by 
NFW, the virial radius r„ of a halo identified at the present 
epoch is related to its virial mass as follows: 



r„ oc M„ 



1/3 



(7) 



From this relation and from eqs. (4), (5), and (6) we obtain 
for the ACDM model that: 



Vm ~ 6.2 X 10" 



M„ 



H-^Mq 



1/3 



Cnfw 
/(cnfwJ 



km s 



(8) 



If Cnfw would not depend on mass, one would have M„ ex 
Vm'^. From the fittings of our halo profiles to the NFW profile 
and using the same definition of virial radius as in NFW, 
we find approximately the dependence of cnfw « M^'^'^''^^ . 
Substituting this dependence in eq. (8) we find that: 



M„ « 5.2 X 10* 



Vm 



km s 



(9) 



This relation is in good agreement with the results obtained 
in our numerical and seminumerical simulations (Fig. 9). 

An intuitive (although only approximate) explanation 
for the halo mass-velocity relation may be given using sim- 
ple scaling relations. For example, Gott & Rees (1975) pre- 
dicted that for objects forming instantaneously (monolithic 
collapse) from density fluctuations with a power-law power 
spectrum the circular velocity and the density of the objects 
should scale as M 12 and M 2 ^ respectively, where n is 
the slope of the power spectrum. For the CDM models at 
galactic scales n « —(2.0 — 2.3). Thus, according to this 
crude analysis the mass scales as Vm'^'^ and Vn'*", respec- 
tively, which is steeper than what we flnd in our simulation. 
The circular velocity in this approximation is the velocity 
at the virial radius of the halo. Note that if the halo den- 
sity does not dependent on the mass, then the mass scales 
as the cube of velocity. In this instantaneous approximation 
the epoch Zc at which the object collapses is related to its 
mass M as (1 -I- Zc) oc M~°', where a = {3 + n)/6, a — 1/6 
for n = —2. However, in the hierarchical formation scenario 
the DM haloes do not form instantaneously; they form in 
a course of aggregation of subunits and accreting material. 
Moreover, due to the random nature of the primordial den- 
sity fluctuations, the MAHs for a given present-day halo 
mass have a dispersion. Nevertheless, one may still define 
a typical epoch of formation of haloes of a given present- 
day mass (Lacey & Cole 1993). For example, this epoch can 
be defined as the average value of the redshifts at which the 
haloes of a given present-day mass M attain half of its mass. 
Using the extended Press-Schechter approximation we cal- 
culate this epoch for several masses in the range of galaxy 
masses and for the ACDM model used here. We obtain: 



l + Zc{M) ocAf" 



(10) 



with a ~ 1/22 — 1/28, i.e. the slope of the collapse redshift- 
mass relation is much flatter than in the case of the in- 
stantaneous collapse (see also AFH98). This implies that 
the densities of the haloes are also less dependent on mass 
than in the case of the instantaneous collapse. Therefore, 
the slope of the mass- velocity relation is smaller than in the 
latter case, and closer to three, which agrees better with our 
numerical results. 

Assuming the spherical top-hat collapse model and as- 
suming that the radius of the virialized object is half of the 
maximum expansion radius, we find that (e.g., Padmanab- 
han 1993): 



Vex {l + Zc^/^M^^^. 



(11) 



Thus, if the mass- velocity relation is of the form M oc V", 
then 



1 + Zc(x M 



2/a-2/3 



(12) 



Comparing this expression with eq.(lO), we find that 
the value for the slope a is ~ 3.2. This is roughly the value 
we obtain in the numerical and seminumerical simulations 
for the ACDM model. In this simplistic analysis we have 
not considered the structure of the haloes; only global scal- 
ing laws were used. Nevertheless, the analysis clearly shows 
that the power-law relation between mass and circular veloc- 
ity (defined at the virial radius) of the haloes is explained by 
the power spectrum of fluctuations and the extended (hier- 
archical) process of formation of the dark haloes. The latter 
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on average also depends on the power spectrum, while the 
scatter in this process (in the MAHs) is determined by the 
statistical nature of the density fluctuation field. 



6.3 Is the Tully-Fisher relation a direct imprint of 
the Mh — Vm relation? 

We address the question of what constraints can be obtained 
by contrasting the observed Tully-Fisher relation (TFR) for 
galaxies and the M — Vm relation for haloes. As it is well 
known, the luminosity in the infrared passbands {H or / for 
example) is a good tracer of the stellar disc mass (e.g.. Pierce 
& TuUy 1992). Therefore, in the assumption that the disc 
mass is proportional to the total halo mass Mh, we would 
expect that the infrared- band TFR is an imprint of the Mh~ 
Vm relation. By comparing the slopes of the Mh — Vm relation 
for isolated halos which we obtained in our numerical and 
seminumerical simulations (§5.2), with the observed TFR 
slopes (e.g., Gavazzi 1993; Peletier & Willner 1993; Strauss 
& WiUick 1995; Willick et al. 1996; Giovanelli et al. 1997), 
we find that indeed this seems to be the case. In other words, 
the Mh/ L ratio in the infrared bands should not depend on 
mass (luminosity). Otherwise the slope of the TFR would 
become different from the slope of the Mh — Vm relation 
(fs 3.2 — 3.3) which already is in good agreement with the 
observational data. 

Regarding the deviations from the Mh — Vm relation, 
they will contribute to the scatter in the TFR. Observational 
estimates indicate a scatter in the TFR of about 0.20-0.45 
magnitudes (e.g., Bernstein et al. 1994; Mathewson & Ford 
1994; Willick et al. 1996; GiovaneUi et al. 1997). Assuming 
again a constant Mh/ L ratio, these estimates correspond to 
a scatter in the Mh - Vm relation of AlogM^ « 0.08 - 0.18. 
These values are in agreement with those we find in our 
numerical and seminumerical simulations. For example, for 
Mh « lO^^/i-iM© we find AlogMh =0.12 and 0.11 for the 
isolated haloes in the N-body and seminumerical simula- 
tions, respectively. As Mo, Mao, & White (1998) and AFH98 
noted, the scatter in the TFR is caused not only by the scat- 
ter in structure of the DM haloes (due to the scatter in the 
MAHs), but also by the dispersion in halo's spin parameter 
A. Nevertheless, Firmani & Avila- Reese (1999b) have shown 
that the quadratic contribution of this latter to the total 
scatter of the TFR is small - only about 25% compared to 
75% contributed to the scatter by differences in the MAHs. 

Our conclusion is that the slope as well as the scatter of 
the Mh — Vm relation of the CDM halos are similar to those 
of the observed TFR and its scatter in the infrared bands. 
This coincidence suggests that the discs formed within the 
CDM halos have a Mh/L ratio in the infrared bands in- 
dependent from mass. There is no room for intermediate 
astrophysical processes (star formation, feedback, gas cool- 
ing) able to introduce a dependence of the infrared Mh/L 
ratio with the mass. Models of galaxy formation and evolu- 
tion where the fraction of the total mass available for form- 
ing stars and the star formation efficiency almost do not 
depend on the total mass of the system, are able to predict 
most of the structural, dynamical and luminous properties of 
disc galaxies, as well as their correlations (Firmani & Avila- 
Reese 1999a,b,c). The observed color-magnitude and color 
TF relations can be well reproduced by these models if the 
luminosity-dependent dust opacity estimated by Wang & 



Heckman (1996) from a large sample of galaxies is intro- 
duced. Wang & Heckman have found that the dust opacity 
of disc galaxies increaeses with their luminosities. This kind 
of correction also might help to match the predicted lumi- 
nosity function in the CDM models with that inferred from 
observations (e.g., Somerville & Primack 1998). 

As a matter of fact, the evolution of the luminous part 
of galaxies is a very complicated process, which goes beyond 
the scope of the present paper. It is obvious that only cool- 
ing gas can produce stars, and, thus, luminosity is defined 
by a complicated interplay between the cooling and heat- 
ing in the baryonic component (e.g.. White & Rees 1978; 
Rees & Ostriker 1977; Blanton et al. 1999; Benson et al. 
1999). But this does not mean that the amount and the 
distribution of the DM are not important. For galaxies with 
1^=100-300 km s~^the total luminosity very likely depends 
on the total mass of the baryons available for star formation. 
The latter correlates with the DM mass. This dependence 
of the luminosity on the DM was observed in hydrodynami- 
cal simulations which include realistic cooling, heating, and 
star- formation processes (e.g., Yepes et al. 1997; Elizondo 
et al. 1999; Steinmetz & Navarro 1999). The correlation ex- 
ists because we are dealing with massive haloes of ym=100- 
300 km s~^for which the gas cools relatively fast on a dy- 
namical time scale and a large fraction of the gas is converted 
into stars. 

The situation is different for haloes with smaller mass 
(Vm < 50 km s^^), which are capable of expelling most of 
their gas if only few supernovae are produced, and which 
can be affected by the intergalactic ionizing background. It 
is also different for larger haloes of Vm > 300 km s~^, which 
host groups or clusters of galaxies. In this case the cooling 
time is long and gas is not converted into stars. Most of 
arguments against a tight L — Vm relation is for those group- 
and cluster- size haloes (e.g., Blanton et al. 1999; Benson et 
al. 1999). In this paper we mostly dealt with galaxy-size 
haloes for which one actually may expect that luminosity 
correlates with Vm- 



7 SUMMARY AND CONCLUSIONS 

We have analyzed the environmental distribution, the outer 
density profiles, and the structural and dynamical correla- 
tions of thousands of galaxy-size DM haloes identified at 
2: = in a cosmological A'^-body simulation of a ACDM 
model. We have also studied and analyzed the formation 
and evolution of DM haloes using an approach based on the 
extended Press-Schechter approximation and on a general- 
ization of the secondary infall model. Our main results and 
conclusions can be summarized as follows. 

1. The density profiles of most of the DM haloes in the 
A-body simulation (typically resolved only down to radii 
0.3 — 0.8rm) are well fitted by the profile given by eq. (2) with 
a distribution of the outer slope l3 such that at the 16%, 50% 
and 84% of the cumulative distribution /3 approximately is 
2.5, 2.9 and 3.9, respectively. The estimated error due to 
small number of particles is large when l3 is large. The slope 
P very weakly anticorrelates with the mass. Our results con- 
firm that the NFW profile shape describes reasonably well 
the intermediate and outer regions of a large fraction of DM 
haloes, particularly the isolated haloes. Our results, how- 
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ever, show that some fraction of haloes have outer profiles, 
which deviate substantially from the NFW shape. 

2. The distribution of the slope j3 and the halo con- 
centration Ci/5 change with the environment in which the 
haloes are embedded. In agreement with previous studies, 
for a given mass we find that haloes in clusters typically 
have steeper outer density profiles and are more concen- 
trated than the isolated haloes. Contrary to a naive expec- 
tation, we find that the haloes in galaxy and group systems 
as well as the haloes with massive companions, systemati- 
cally have fiatter and less concentrated density profiles than 
the isolated haloes. The fact that the halo density profiles do 
not follow a continuous trend along the cluster-group-field 
sequence suggests that the difference between clusters and 
groups cannot be viewed only as a question of density. 

3. Approximately 70% of the galaxy-size DM haloes of 
130 km s"'^ < Vm < 350 km s~^ are very isolated systems 
in the sense that they are not contained within larger haloes 
and they do not have massive companions {V^'^'^ > 0.7Vm 
or M™™'' > 0.3M/i) within a radius equal to 3 times their 
own radii. The « 13% of the haloes are contained within 
larger haloes. The haloes in pairs or multiple systems con- 
stitute the « 17% of all haloes. 

4. The parameter Ci/5 is a good estimator of the halo 
concentration, independent of the profile fitting. The less 
massive haloes tend to have larger values of Ci/5 than the 
more massive haloes. 

5. The galaxy-size haloes exhibit a relation between 
their masses and maximum circular velocities, Mh oc V^, 
with a ~ 3.3 and a ~ 3.5 for the isolated and cluster 
haloes, respectively. This relation may be considered as an 



I 



imprinu of the primordial density fiuctuation field. For a 



mass ot W^^ h~^ Mq the rms fractional velocity deviation 
cry/ < Vm > from this relation is ~ 0.085 for isolated haloes 
and ~ 0.128 for cluster haloes. The deviations correspond to 
logAMj, ~ 0.11 and ~ 0.18 for isolated and cluster haloes 
respectively. The deviations of the Mh — Vm and Ci/5 — Mh 
relations are tightly correlated. For a given mass the more 
concentrated haloes have larger Vm- 

6. The distribution of the parameter (3 obtained with 
the seminumerical approach, is similar to the distribution of 
/3 for isolated haloes in the A'^-body simulations; the median 
is at /3 ~ 2.78. The Mh — Vm and Ci^^—Mh relations and their 
dispersions are similar to those of the isolated haloes, too. 
This affrppmeut between two completely different methods 



IS enco mraging 



To conclude, we have shown that the shapes and con- 
centrations of DM haloes exhibit a diversity and systematic 
dependence on the halo's environment, the NFW shape be- 



ing close to the the average shape and concentration. The di- 
versity mid dependence on eiiviiuniueiiL can be iiupuiLaiiL in 
shaping the properties of galaxies and their scatter. There- 
fore, studies of galaxy formation and evolution should make 
an atte mpt to account for these effects. 
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Table 1 . Enviromental distribution of DM haloes with maximum 
circular velocities Vm > 130 km s^^. 



Environment 



Number of haloes Percentage of haloes 



elongs to: 




Cluster 


227 


Group 


112 


Galaxy 


98 


ot in a larger halo: 




Isolated 


2456 


Pairs 


605 



Total 



3498 



6.5 
3.2 
2.8 



70.2 
17.3 

100 
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Figure 1. Differential and cumulative distributions of the dis- 
tance from the center of distinct halo (not contained within 
a larger halo) to its nearest significant companion {V^™^ > 
O.TVm)- The distance is scaled to the radius rj, of the halo. The 
bin width is Admin = 0.67. 



Figure 3. Differential distribution of the slope /3 for all the haloes 
with 130 km s^^ < Vm < 350 km s^^ (solid line) and only for 
the haloes with the low values of x^ of the fit (thin solid line). 
The arrows indicate the 16%, 50%, and 84% of the cumulative 
distribution of all the haloes. The bars show the errors of /3. 
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Figure 2. Fraction of isolated haloes in the halo catalog as a func- 
tion of the lower limit on the circular velocity of the companion 
fy = V^^'^ /Vm- The isolated haloes are haloes not contained 
within larger haloes and without any companion with circular 
velocity V^"^^ > /yVm within 3ri^, where Vmand r/j are the 
circular velocity and radius of isolated halo. The lines are for dif- 
ferent samples. The only difference between these samples is the 
lower limit on Vm (shown in the panel) allowed for the isolated 
haloes. The curves are truncated at the value of fy at which the 
sample becomes incomplete in the sense that the companions of 
the smallest isolated haloes can be smaller than 90 km s^^. 



Figure 4. Same as Fig. 3 but for different environments indicated 
in the panels, and for the haloes obtained in the seminumerical 
simulations (dashed line in the upper panel). For the isolated 
haloes (upper panel) the distribution for haloes larger than the 
galaxy sizes, is also shown (dotted line). The bin widths A/3 used 
to calculate the distributions shown in the upper, medium, and 
lower panels were fixed to 0.20, 0.38 and 0.30, respectively. The 
bars plotted for the distribution of the group and cluster haloes 
show the error we estimate in the determination of the param- 
eter /3 when /3 ai 2.5 and 3.9 for the group and cluster haloes, 
respectively. 
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Figure 5. Density profiles of lialoes in different environments 
(crosses) and tlie fittings to tliese profiles using the generalized 
NFW profile given by eq. (2) (solid lines). For each sample three 
haloes were randomly chosen in three mass ranges with /3 around 
the peak of the distribution of the corresponding sample (/3 is 
shown in each panel). For the haloes in clusters we have chosen 
/3 = 4 instead of f3 at the peak of the distribution. The stright 
lines indicate different slopes. 



Figure 7. The concentration parameter c^m vs. the concentra- 
tion parameter c = Tf^/rs, where r^ is the scale radius in the 
fitting formula eq. (2), for all the haloes with /3 Ri 2.5 (stars), 
fs 3.0 (dots), and Ri 4.0 (crosses). 
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Figure 6. The goodncss-of-the-fit parameter •v/x'^/^bins 
centages for the cases when the halo density profiles were fitted 
to the Hernquist and NFW profiles. The upper and lower panels 
are for the isolated haloes and the haloes in clusters, respectively. 



Figure 8. Dependence of concentration c^/s on halo mass M^ 
for haloes in different environments (indicated in the panels), and 
for the haloes from the seminumerical simulations (panel d) . The 
solid lines refer to the average concentration calculated for several 
mass bins. The logarithmic widths of the bins are AM;, = 0.10 — 
0.12. The standard deviations are represented with the dashed 
lines. The thin solid lines are the linear regressions applied to all 
the haloes of the given sample. 
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Figure 9. Dependence of mass Af^ on the maximum circular 
velocity Vm for the isolated haloes (a), the haloes in clusters 
(b), and the haloes obtained in the seminumerical simulations 
(c). The same line code of Fig. 8 is used. The slopes of the linear 
regressions (thin solid lines) are indicated within each panel. 



Figure 11. Correlation among the residuals of the M^ — Vm and 
Cj /5 — Mf^ relations for the isolated haloes (upper panel) and for 
the haloes in clusters (lower panel). 
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Figure 10. Fractional rms scatter in velocity of the mass- velocity 
relation as a function of the mass. The solid and short-dashed lines 
are for the isolated and cluster haloes, while the long-dashed line 
is for the haloes in the seminumerical simulation. 
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